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1. Introduction 

Scattering phase shift is an important physical observable in its own right. However, for nu- 
clear physics, it is more advantageous to convert it into a form of the nuclear potential. Once such 
a nuclear potential is at our disposal, we can conveniently use it to study a variety of nuclear phe- 
nomena based on the effective degrees of freedom, i.e., the nucleons. It provides us with physics 
insights into the structures and the reactions of atomic nuclei, as well as the supernova explosion 
of type II and the structure of neutron stars through the equation of states of cold and dense nuclear 
matter. Enormous efforts in this line are integrated into a form of realistic nuclear potentials By 
using about 40 — 50 adjustable parameters, they can reproduce several thousands of experimental 
NN data with ;^^/NDF ~ 1, which consist of the scattering phase shifts and the deuteron property. 
Also, the potentials from the chiral effective field theory attract an growing interest 

Unlike the nucleon sector, only a limited number of experimental information is available in 
the hyperon sector. This is due to the absence of accelerator facilities, which can generate direct 
hyperon beam. If there were a sufficient number of experimental information on the hyperon- 
hyperon and hyperon-nucleon scatterings, they would be converted into realistic hyperon potentials, 
which could help us study the structure of hyper-nuclei and possible generation of hyperon matter 
in the neutron star core. 

The standard method to obtain the scattering phase shifts in lattice QCD is Luscher's finite 
volume method [^]. It can be used to provide QCD predictions/postdictions for the scattering phase 
shifts not only in the nucleon sector [Q, §], but also in the hyperon sectors One may come up 
with a straightforward way to obtain realistic inter-baryon potentials in lattice QCD, i.e., sufficient 
number of scattering phase shifts are generated by Luscher's method at the initial stage, which are 
converted to inter-baryon potentials with the help of the inverse scattering theory. However, this 
is difficult in practice, because it involves an infinite number of scattering phase shift at the initial 
stage. Thus, it is desirable to have a direct method to obtain realistic inter-baryon potentials in 
lattice QCD. 

The method recently proposed by Ref. is such a method. It can be given a background 
in terms of Luscher's finite volume method. By using the effective Schrodinger equation, it con- 
structs nuclear potentials from the Bethe-Salpeter (BS) wave functions generated by lattice QCD. 
Since the information of the scattering phase shift is embedded in the long distance part of the BS 
wave functions, it is possible to generate a realistic nuclear potential, which reproduces the QCD 
predictions of the scattering phase shift extracted by Luscher's method. 

In this paper, after a brief review of the general idea how to construct a realistic nuclear po- 
tential in lattice QCD, we present quenched QCD results of the central and the tensor potentials 
obtained at the leading order of the derivative expansion. After the discussion of the convergence of 
the derivative expansion, we give dynamical QCD results by using 2-1-1 flavor gauge configurations 
generated by PACS-CS Collaboration. Finally, we mention the hyperon potentials {N'E and A'^A), 
to which our method can be equally applied. 

2. General idea to construct nuclear potential in QCD 

We consider (equial-time) Bethe-Salpeter (BS) wave function for two nucleons in the center 



2 



Lattice study of nuclear forces 



Noriyoshi Ishii 



of mass frame. By choosing particular composite interpolating fields p{x) and n{x) for proton and 
neutron, the BS wave function is defined as 

Wki^-y) = {^W)n{y)\pCk)n{-k)M) = }^^{0\tW,^^^^ (2-1) 

where k denotes the asymptotic momentum of the proton relative to the neutron. It is related to 



the relativistic total energy of the state as Pq = 2y m^ + k^ with m^? being the nucleon mass. Quite 
naively, this matrix element may be regarded as an amplitude to find three quarks at x and another 
three quarks at y, where p{x) and n{y) are used to probe nucleons in the state \p(k)n{— k) , in) . 
At long distance |jc — j| — ?■ large, it shows an asymptotic behavior, which is characterized by the 
scattering phase shift 5/ {k) in exactly the same way as a scattering wave function in the quantum 
mechanics as ||8|, ^ [Io| ] 

?,m{k\x-y\-nl/2 + 8i{k)) 
k\x — y\ 



^r^^{x-y)^A ^"^^'^'■^ ^IJ-'i;^''^^''" +.... (2.2) 



To prove this behavior, Nishijima-Zimmerman-Haag (NZH) reduction formula [ 1 1 1 is conveniently 



used. Note that any local composite nucleon field N{x) leads to the same asymptotic behavior 



Eq. ( |2.2[ ), as far as N{x) has non-vanishing overlap with a single nucleon state, i.e., (0|A'^(x)|A'^) ^ 0. 
In Luscher's finite volume method, which is the standard method to calculate the scattering phase 
shift in lattice QCD, the phase shift embedded in BS wave function in this manner is extracted from 
the energy spectrum in a finite periodic box (For explicit use of BS wave function in Luscher's 
finite volume method, see Ref. [^.) 

For nuclear physics, it is more advantageous to convert the data of the phase shift into a form 
of nuclear potentials. We therefore wish to extend LUscher's method so as to obtain the nuclear 
potentials directly. For this purpose, we use the remarkable similarity in the asymptotic behaviors 
between the BS wave function Eq. ( p^ ) of QCD and the scattering wave function in the quantum 
mechanics. This similarity motivates us to construct a nuclear potential so that it can reproduce all 
the BS wave functions simultaneously in wide range of energy region. Then, the resulting potential 
can reproduce the phase shifts predicted by QCD. In this way, it becomes possible to construct 
realistic nuclear potentials by lattice QCD. 

To proceed, we define the nuclear potential U {7,7') by the effective Schrodinger equation [ IC] 

{A + k^)Yi{7)=mM J d'r'U{7,r)Yi{7'), (2.3) 
where m^r denotes the nucleon mass. We demand this equation to be simultaneously satisfied by 



1//^ (7) in wide k region (or wide energy 2 y mjj + k'^ region). Note that U (7, 7') is most generally a 
non-local potential, and that, with our definition, U if^f) does not depend on the relativistic total 



energy Po = 2y + ^2. 

Several comments are in order. 

(i) With our prescription, precise forms of potentials depend on particular choices of interpo- 
lating fields p{x) and n{x). However, even if their particular shapes are different, these potentials 
lead to the same phase shift. Remember that these potentials are constructed so as to give the phase 
shift obtained by Liischer's method, which does not depend on a particular choice of interpolat- 
ing fields. The situation is analogous to the unitary transformation in quantum mechanics, i.e.. 
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shapes of potentials are changed by unitary transformations without affecting any observables. The 
potentials, which lead to the same phase shift, are said to be phase-shift equivalent to each other. 

(ii) Here, for simplicity, we do not pay attention to the orthogonality of the BS wave functions. 
In general, the BS wave functions are not expected to be orthogonal to each other with respect 
to the three-dimensional integral as j d^x VA|,^^(^)i/A^.^^(jc) / ^ 8^ (k — k'). If the violation of 
the orthogonality is serious, we have to take care of the orthogonalization to avoid non-hermitian 



potentials [10] 



3. The Derivative Expansion 



To construct the non-local potential U{r,r') in Eq. ( |2.3| ), it is necessary to generate infinite 
number of BS wave functions, which is difficult in lattice QCD. The reason is two fold, (i) The 
energy spectrum is discretized in a finite spatial box. (ii) Although there is a method to access 
excited states in lattice QCD, it becomes the more difficult to access the higher excited states. 
Thus, we need an approximation, which enables us to construct U{x,^) with a limited number 
of BS wave functions. For this purpose, we use the derivative expansion. We can start with the 
leading local potentials, and then take into account possible non-local terms (potentials, which 
contain derivatives) order by order. If the non-locality effect appears to be large, the convergence 
can be improved by changing the interpolating fields p{x) and n{y) in the sink side. 

To proceed, we impose general requirements on the non-local potential ?7(x,x') arising from 
the translational invariance, Galilean invariance, symmetry condition (identical particle condition), 
spatial rotation, spatial reflection, time-reversal invariance and hermiticity. The most general form 
has been derived in Ref. JT2|], to which we apply the derivative expansion. We are left with 



U'{xX) = V{,^{xy)5\x-^) (3.1) 
V^^(f,V) = vl{r) + vUr) Oi ■ a2 + Vi{r) Sn + V[s{r)L-S + OiV^), 

where / indicates the total iso-spin of the two nucleon system. Oi and 02 act on the spin indices of 
the first and the second nucleons, respectively. Sn = 3(ai ■x){d2 ■x)/x^ — Oi - 02 is referred to as 
the tensor operator. L = ixxV denotes the orbital angular momentum operator. 5 = (ai + a2)/2 
denotes the total spin operator. Since d\ ■ 02 reduces to 1 for spin triplet, and —3 for spin singlet, 
Vo{r) + V(j{r) di • 02 is conveniently combined into the form of the "central potentiaF Vc{r) as 

Vc{r,'So)=Voir)-3Vair), Vc{r,^ S^) = Vo{r) + Va{r). (3.2) 

Vrir) and Vz,s(r) are referred to as the tensor potential and the LS potential, respectively. Note that 
Vc{r), VV(r) and VLs{r) play important roles in conventional nuclear physics. 

At the leading order, we truncate the nuclear potential by neglecting all the derivative terms as 

VNNix,y)=Vcir)+VTir)Sn + 0{y). (3.3) 



We insert this into the effective Schrodinger equation Eq. (2.3) to have 



'— + Vc{r) + VT{r)Sn) ¥%{7) = — V^j(r). (3.4) 
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Note that the reduced mass of two nucleon system is /i = mf^jl. For = 0+ 1=1 channel, 

a further reduction is possible. Since acts as zero on wave functions in channel, we can 
rewrite it as 

Vc{r; Sq) = — + f—. (3.5) 

We use this formula to obtain the central potential in channel at the leading order of the deriva- 
tive expansion. For = 1+ (^Si —^Di coupled system) 7 = channel (the deuteron channel), the 
procedure is slightly involved, which is to be considered in Sect. ^. 

At the next to leading order, we include the terms, which contains a single derivative, i.e., 
VLs{r)L-Sas 

Vnn{x,V) = Vc{r)+VT{r) Sn + VLs{r)L-S + 0{V^). (3.6) 



This is inserted into the effective Schrodinger equation Eq. ( |2.3| ). Note that the action of L • 5 on 
= 0+ ('^o) channel vanishes. Therefore, the formula Eq. (3.5) does not change at this order. 



In contrast, L- 5 gives a non-vanishing contribution to J = 1+ ( — Di) channel. Hence, the 
formula to calculate Vc(r) and Vjir) is modified. We need an additional BS wave function to obtain 
these three potentials at this order. 

At the next to next to leading order, we include O(V^) terms in the potential, which is inserted 



to the effective Schrodinger equation Eq. (2.3). To obtain these potentials, we need further BS 
wave functions. We perform this procedure repeatedly to obtain higher derivative terms by using 
increasing number of BS wave functions. 

It is important to examine the convergence of the derivative expansion. The non-local potential 
is faithful to the scattering data in wide range of energy region, while it may not be so after the 
derivative expansion is applied. If the convergence appears to be unsatisfactory, improvement has 
to be done by changing interpolating field of nucleon. 



4. Lattice QCD setup 

We use quenched QCD unless otherwise indicated. We employ the standard plaquette gauge 
action with j8 = 5.7 to generate gauge configurations on the lattice of the size 32^ x Nt (Nt = 32 and 
48). The scale unit is introduced by rho meson mass in the chiral limit, which leads to the lattice 
spacing a^^ = 1.44(2) GeV (a ~ 0.137 fm) [Q]. The spatial extension amounts to L = 32a ~4.4 fm. 
1000-4000 gauge configurations are used in our calculations. Quark propagators are generated by 
employing the standard Wilson quark action with the hopping parameters K = 0.1640,0.1665 and 
0.1678, which correspond to ~ 731,529 and 380 MeV and mn ~ 1558, 1334 and 1197 MeV, 
respectively. Unless otherwise indicated, Dirichlet and periodic boundary conditions are imposed 
on quark fields along the temporal and spatial directions, respectively. 

To obtain the BS wave function, we generate the four point correlator of the nucleon field as 

G{x-y,t) = {0\T [p{x,t)n{y,t)p{t = 0)n{t = 0)] 1 0> , (4.1) 

where p{x) and n{x) denote local composite fields for proton and neutron, for which we employ 
the standard ones as p{x) = Satc {u^{^)Cy5dh{x)) u{x), n{x) = eatc {'^l {^)CY5db{x)) d{x). These 
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fields are represented by Heisenberg picture in imaginary time. p(f) and n{t) denote interpolat- 
ing fields for proton and neutron (wall source), i.e., p(f) = Y.x.y.z^abcU-c{x.,t) {db{y j)CysUa{Z:tY) > 
n{t) = Y,x,y,z^abcdc{^j) {dbiy,t)CY5Ua{z,tY) ■ Notc that the total spatial momentum of the system 



vanishes because of the wall source. Therefore Eq. ( [4. 1| ) becomes a function of x — y due to the 
translational invariance. Eq. ( p~T| ) for the large Euclidean time f > is dominated by the contribu- 
tion from the lowest-lying state as 

G{x-y,t) = £ (0 \p{x,t)n{y,t) \ E,,) (£„ |p(0)n(0)| 0> = Aq e'^^' WE,{x-y) + • • • , (4.2) 

n 

where £0 denotes the energy of the lowest lying state j^o), and Aq = (£'o|p(0)n(0)|0). ^/^^(x — y) = 
(0\p{x)n{y)\EQ^ denotes the BS wave function. Needless to say, YEo{r) does not depend on a 
particular choice of the interpolating fields in the source side, as long as (SqIp n|0) 7^ 0. (We use 
the wall source for the efficiency reason.) Quantum numbers of the wave function such as are 
controlled by quantum numbers of the interpolating fields in the source side. Because we adopted 
the wall source, = 0+ and 1+ are obtained by combining the spins of p and n. 

5. Central potential at the leading order of the derivative expansion 
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Figure 1: (a) BS wave functions in ^Sq and ^Si channels for m^: ~ 529 MeV. The inset is a 3D plot of 
y/{x,y,z — 0) for ^Sq. The curves denote the results of the fits using the Green's function of Helmholtz 
equation in the region lla<r< 16a and (b) The central potential in channel and the effective central 
potential in ^Si channel for ~ 529 MeV. 

Fig. [I](a) shows the quenched result of BS wave functions in ^Sq and ^Si channels. To pick 
up "s-wave" component (^5i) from the BS wave function in 7^ = 1+ channel, which is ^Si-^Di 
coupled system, we make a spatial average with respect to the cubic group O as Yapi^'^'^ Si) = 
^T^geoVajiis^)- Calculations are fully performed for r < 0.7 fm, while, for r > 0.7 fm, we 
restrict ourselves to the points on the coordinate axes and their nearest neighbors to reduce the 
calculational cost. (A rapid change of the potential is expected for r < 0.7 fm, whereas a rather 
mild change is expected for r > 0.7 fm.) We see that there are shrinks at short distance, which 
suggest the existence of repulsion. By using Eq. (3.5), the central potential Vc{r) in ^Sq channel is 
constructed at the leading order. The result is shown in Fig. |l]^b). 
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To obtain the central potential in channel, it is necessary to consider a coupled Schrodinger 
equations of and channels, which will be discussed in Sect. ^. Here, we simply apply the 
same formula Eq. (3^) to the wave function in ^Si channel. Note that the resulting central potential 
can reproduce the wave function without involving the tensor potential. Such a central potential 
is referred as the "effective central potential", in which the effect of the tensor potential is embedded 
implicitly. The result is also shown in Fig. [I](b). 

In Fig. |l|(b), we wee that phenomenological properties of the central nuclear potentials are 
reproduced. A repulsive core at short distance is surrounded by an attraction at medium distance. 
The effective central potential in channel tends to be more attractive than the central potential 
in '^o channel. This is desirable for the existence of a bound state (deuteron) in channel in 
reality. (No bound state exists in '^o channel.) 

The non-relativistic energy E = Ir/niiM in Eq. (3^) is obtained by making a fit of BS wave 



functions with Green's function of Helmholtz equation defined as 

{A + k^)G{x;k^) = -5L{x), 



(5.1) 



where 5l{x) = Y^nei? ^^P {^Ttin-x/L) denotes the periodic delta function in the three dimensional 
torus of spatial extension L. The fit is performed in the region 1 la < r < \6a, where the interaction 
is seen to become negligible from a plot of Ai/a(.x) / ^f{x) |^. The resulting non-relativistic energies 
are quite small, i.e., E = -0.509(94) MeV for channel and E = -0.560(110) MeV for 
channel. Note that negative E does not necessarily mean a formation of a bound state. This is 
because two nucleons cannot be separated from each other beyond the range of interaction in a 
finite volume. The result is also shown in Fig. |I|(a), where Green's functions along the coordinate 
axis are plotted with solid lines. 
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Figure 2: (a) The central potential in channel for three quark masses and (b) those with r^ multiplied. 

Fig. ^a) shows the central potentials in 'So channel at the leading order for three values of the 
quark mass. To see the strength of the potentials more accurately, we plot them with a factor in 
Fig. ^b) to reflect the effect of three-dimensional volume element. We see that, as the quark mass 
decreases, the repulsive core at short distance grows rapidly, and the attraction at medium distance 
is enhanced gradually. The enhancement of the attraction is natural, since the smaller quark mass 
enables the virtual pion to propagate the longer distance. These results suggest the importance of 
direct lattice QCD calculations in the light quark mass region. 
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6. Tensor Potential at the leading order of the derivative expansion 



In order to construct the central and the tensor potentials, it is necessary to take into account the 
coupling of s-wave(Z = 0) and d-wave(Z = 2) components of BS wave function in 7^ = 1+ (/ = 0) 
channel. Note that the coupling is induced by the tensor potential, due to which the deuteron is 
generated as a bound state in reality. The tensor potential plays a key role in the stability of atomic 
nuclei and the saturation of nuclear matter density. Its importance at short distance is pointed 
out recently by experimental studies of Short Ranged Correlated (SRC) nucleon pair, which may 
affects the structure of the cold dense nuclear system such as neutron stars [|l^]. However, the 
experimental determination of tensor potential at short distance is difficult, because it appears on 
top of the repulsive core and the centrifugal barrier. 

Construction of tensor potential depends on the quality of d-wave wave function on the lattice. 
On the lattice, we construct a BS wave function in representation of the cubic group, which 
corresponds to 7^ = 1+ up to 7 > 4 contamination. We decompose it into orbi tally A\ part (/^('^^(r), 
which corresponds to s-wave up to / > 4 contamination, and orbitally non-Aj^ part i/A(^)(r). The 
decomposition is performed by using the projection operators P and Q defined as 



(6.1) 



where O denotes the cubic group with 24 elements. From the orthogonality relations of the charac- 
ters of the representations, we find that the orbital part of xf/^^^r) consists of either or 
representations. £"+ and correspond to d-wave up to / > 4 contamination. Orbital represen- 
tation corresponds to g-wave (/ = 4) up to / > 6, which enters through = 4+ component of 
through the relation 4+(7^) = l(spin) (8)4+ (orbital). Fig. |3|(b) shows BS wave functions for T^^ 
(J^ ~ 1+) and the azimuthal quantum number M = 0. We see that i/A^"^) is single-valued, whereas 
y/^(^)'s are multi-valued. Since the angular dependence manifests itself as multi-valuedness, it 
follows that Y^^^ is dominated by s-wave contribution. To consider \lf^^\ we note that d-wave 
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Figure 3: (a) BS wave functions on the lattice and (b) BS wave function after removing the spinor harmonics 
factors. 
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component is proportional to the ''spinor harmonics" in S0(3) symmetric limit as 

2 



VAio(r) 1/^11 (r) 



y('-2) (t.\ 




2 

'^/6 



J'2,o(r) 



y2.+i{7) 



(6.2) 



To examine whether \\r^^^ is dominated by d-wave or not, we divide i//^'^' by these spinor harmonics 
factors. Fig. ^(b) shows the results. We see that i//^'^' becomes single valued, which indicates that 
is dominated by d-wave. 

To separate the s-wave and the d-wave parts, we apply the projection operators P and Q to 



Eq. (3.4). Since Vcir) and VTir) commute with P and Q due to the rotational invariance, we have 

A 



A 



Pn{r)+Vc{7)Pnir)+VT{r)PSu¥%{r) = —P¥% 
QYkir)+Vc{m¥i{r)+VT{mSn¥k(r) = —Qni^)- 



(6.3) 



Note that each of these two equations has two spinor indices. For upper line, we have essentially a 
unique choice. ((0, 1) and (1,0) components agree to each other up to an overall sign.) In contrast, 
we can play with a particular choice of spin components for the lower line. Eq. ( |6!2| ) suggests that 
(0,1) and (1,0) component correspond to £ -representation, whereas (0,0) and (1,1) component 
correspond to r2-representation. Since S0(3) symmetry is not exact, the results depend on how 
we choose d-wave wave function. For simplicity in this section, we choose (0, 1) spin component 
from the lower line for d-wave. Eq. (|6.3[) is arranged as 



Vcir) - — 

J niN 

Vjir) 



(6.4) 



which can be algebraically solved for Vc{r) and Vrir) point by point. Unlike the central potential, 
the tensor potential does not involve an additional shift by — which adjusts zero at the spatial 
infinity. In Sect. Isl we employ another choice for the d-wave from the lower line by combining the 



four spinor components with the spinor harmonics Y(r) = Y 



(1=2) 
/=1,M=0 



(r) as 



(6.5) 

Once S0(3) is realized as a good symmetry, this becomes the best choice. 

Fig. Q^a) shows the results of the tensor potential together with the central and the effective 
central potentials. The shape of our tensor potential is similar to the one-boson exchange result, 
which is obtained by the cancellation between the pion exchange and the rho meson exchange [0]. 
The difference between the central and the effective central potentials is understood by treating the 
tensor potential in the second order perturbation theory. It seems to be smaller than phenomeno- 
logically expected. This is due to the heavy quark mass. Indeed, as is seen in Fig. ^b), the tensor 
potential is enhanced in the light quark mass region. 

A spike in the tensor potential at r ~ 0.5 fm is due to the zero of F2.o(^) °^ 3cos^ 6 — \ . Note 
that i2,o('^) vanishes on the lines r= (in, in, ±.n). In the vicinity of these lines, it becomes difficult 
to solve the coupled equation numerically, which leads to the accumulated statistical error. 
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Figure 4: (a) The tensor, the central and the effective central potentials in ^Si — ^ Di coupled channel for 
nipi = 529 MeV and (b) the tensor potentials for three values of quark masses. 



7. Convergence of the derivative expansion 

We examine the convergence of derivative expansion by comparing two potentials generated 



at two energies £ ~ MeV and £ ~ 45 MeV []15|]. As soon as the derivative expansion is truncated, 



energy dependence may appear in the potential, because it is not in general possible to reproduce 
all the BS wave functions simultaneously only with truncated degrees of freedom. In this way, the 
energy dependence of the potential is related to the non-locality of the potential, which makes it 
possible to check the convergence of derivative expansion by examining the energy dependence. 

We generate two potentials by imposing different spatial boundary conditions on quark fields. 
A potential at £" ~ MeV is generated with the periodic boundary condition (PBC). A potential at 
£■ ~ 45 MeV is generated with the anti-periodic boundary condition (APBC). With APBC, since 
each nucleon consists of odd number of quarks, a nucleon is also subject to APBC, so that its 
spatial momentum is discretized as pi = {Ini + l)7i/L with n,- G Z. For a two nucleon system, the 
interaction (nuclear force) induces a modification from its free value, i.e., p^'''^'^ ~ (2?i,- + \)n/L. 
Note that the smallest spatial momentum for APBC is /J^''^') ~ (ib7r/L,ib7r/L,ib7r/L), which does 
not vanish. In the box with L ~ 4.4 fm, |^(''''''| ~ \/37r/L = 244 MeV. Comparison is made by 
using the setup with niji = 529 MeV, m/v ~ 1333 MeV, which leads to the non-relativistic energy 
of the lowest-lying state E = k^/mN ^ 45 MeV. The results are shown in Fig. ^. We see that the 
agreement is quite good except small deviations at short distance. The structures appearing in the 
region r > 1 fm for APBC turn out to be caused by a small contamination of an excited state [|l^]. It 
follows that the derivative expansion works, and that our local potential constructed at the leading 



order can be reliably used in the energy region £ ~ — 45 MeV (See Ref. [15] for detail.) 



8. 2+1 flavor QCD result of nuclear force with PACS-CS gauge configuration 

In order to study the quantitative features of nuclear potentials, it is necessary to resort to 
dynamical QCD performed in the light quark mass region employing a large spatial volume. 
PACS-CS Collaboration is generating such gauge configurations, i.e., 2-1-1 flavor gauge config- 
urations, which cover the physical quark mass employing large spatial volumes L ~ 3 — 6 fm 
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Figure 5: Comparisons of leading order potentials generated at two energies, i.e., £ ~ MeV by PBC and at 
£ ~ 45 MeV by APBC. The left, middle, right figures show the central potential in '^o channel, the central 
and the tensor potentials in —^D\ coupled channel, respectively. 



[16]. We use PACS-CS gauge configurations to obtain 2+1 flavor QCD results of nuclear poten- 
tials. The gauge configurations are generated by employing Iwasaki gauge action at j3 = 1.90 
on 32^ X 64 lattice and 0(a)-improved Wilson quark (clover) action with a non-perturbatively 
improved coefficient csw = 1.715. nij^, mk and ma are used to determine the scale unit a^^ = 
2.176(31) GeV (a ~ 0.091 fm) leading to the spatial extension L = 32a ~ 2.90 fm 0]. To cal- 
culate nuclear potentials, we use three series of PACS-CS gauge configurations with {K^d, K'i) = 
(0.16700,0.16400), (0.16727,0.16400) and (0.16754,0.16400), which correspond to m„ ~ 701, 
570 and 41 1 MeV and mn ~ 1583, 1412 and 1215 MeV, respectively. 
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Figure 6: (a) 2+1 flavor QCD results of the central and tensor potentials for = 101 MeV, and (b) 
quenched QCD results of the central and tensor potentials for nij^ = 731 MeV. 

Fig. I^a) shows the 2+1 flavor QCD results of the nuclear potentials for m;;: ~ 701 MeV, which 
should be compared with the quenched results in Fig. ^(b) of comparable pion mass niji ~ 731 
MeV. We see that the repulsive cores at short distance and the tensor potential become significantly 
enhanced. The attraction at medium distance tends to be shifted to outer region, whereas it remains 
almost unchanged in magnitude. Although these changes may be caused by dynamical quarks, they 
may be due to a lattice discretization artifact. We need further information to conclude. 

Fig. ^ shows the central potential Vc{r) in channel and the central potential Vc{r) and the 
tensor potential VV(r) in —^Di coupled channel, together with those multiplied. Similar 
tendencies are observed as the quenched QCD, such as the enhancements of the repulsive cores at 
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Figure 7: 2 + 1 flavor QCD results of nuclear potentials for three values of quark mass, (a) Vc{r) in '5o 
channel, (b) Vc{r) in D\ coupled channel, (c) V-^ir) for —^D\ coupled channel, (d) r^Vc{r) in '^o 

channel, (e) P'Vcir) in -3 D\ coupled channel, (f) r^VT(r) for ^Si —^D\ coupled channel. 



short distance, the attractive pockets at medium distance, and the strength of the tensor potential. 




Figure 8: (a) Scattering phase shifts in So channel from the lattice nuclear potentials and (b) those with the 
empirical one. 

Fig. ^ shows the phase shifts obtained from the nuclear potentials (not by Liischer's method). 
At low energy, the phase shift grows up, which is caused by the attraction at medium distance. At 
high energy, the phase shift decrease, which is a consequence of the repulsive core at short distance. 
One may wonder why the order of the phase shifts of m^i = 570 and 701 MeV is inverted. The 
reason seems to be that the repulsive core grows more rapid than the attraction grows. Qualitative 
shape of the phase shift is seen to be reasonable, which is due to the fact that the qualitative features 
of the nuclear potential are already reproduced. However, comparing with the empirical one, the 
strength is not satisfactory at all, which suggests the importance of the light quark mass effect. 

A technical comment is in order. To obtain the scattering length, the ground state saturation 
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has to be achieved to an accuracy of, at most, around 1 MeV, which is about 0.05 % of the total 
mass of the two nucleon system. For the moment, such a high precision is not yet attained in our 
calculation. Significantly large t and, accordingly, the large statitics are required. To calculate 
only the scattering length, the smeared source is better than the wall source for such high precision 
calculations. The smeared source, with the help of the average over the relative coordinate x in 
the sink side, projects out the excited state contamination from the temporal correlation. However, 
to calculate the BS wave function, which measures the spatial correlation, the sink has to be un- 
averaged, and the ground state saturation has to be achieved point by point uniformly in the spatial 
directions. Note that the smeared source creates a spatially squeezed BS wave function in the small 
t region, which gradually broadens during the temporal evolution, until the ground state shape is 
achieved. Here, the convergence in the region \x\ > 1 fm is quite slow, and unreasonably large t 
is required for the uniform saturation by the ground state. For future applications to the nuclear 
physics, it is necessary to seek for a better source, which makes it possible to achieve the uniform 
saturation of the ground state BS wave function more efficiently. 



9. Hyperon potentials 



Hyperon potentials (hyperon-nucleon and hyperon-hyperon) serve as the starting point in 
studying the hyper-nuclei structure. They have large influence on the hyperon matter generation in 
neutron star core. In spite of their importance, we have only a limited knowledge of the hyperon 
potentials, because of the lack of experimental information. Since we do not need any information 
from the scattering experiment, we apply our method to construct the hyperon potentials. The first 
attempts have been made to construct A'^S (/ = 0) potentials and A^A. The results are shown in Fig. ^ 



(See Refs. [ ]17| , |18| ] for detail.) These calculations are being extended to Nl. and AA potentials. It 
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Figure 9: NZ potentials(/ = 1) from quenched QCD (left) and A^A potentials from 2+1 flavor QCD (right). 

is important to examine the convergence of the derivative expansion. After this process, the reli- 
ability of the potential is guaranteed in the well-defined applicability region. Such potentials may 
be used as alternatives to experimentally constructed potentials, before experimental information 
on various hyperon scatterings becomes fully available. 

Hyperon potentials give us another benefit. Their flavor structure is expected to provide us 
with an important key to unveil the physical origin of the repulsive core. To obtain a simplified 
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picture, investigations in the flavor SU(3) idealized limit would be useful and interesting. Works 
are in progress along this line [[l9|]. It would be also interesting to compare these numerical results 
with analytical ones obtained by the operator product expansion [20]. 



10. Summary 

We have reported lattice study of nuclear potentials based on the equal-time Bethe-Salpeter 
(BS) wave function for two nucleon system. We have defined the nuclear potential by the effective 
Schrodinger equation demanding that it should generate the BS wave functions in wide range of 
energy region simultaneously. The method generates realistic nuclear potentials, because of the 
remarkable similarity in the asymptotic behaviors between the BS wave function of QCD and 
the non-relativistic wave function of scattering state in the quantum mechanics. The central and 
the tensor potentials have been obtained at the leading order of the derivative expansion, which 
show the qualitative features of the phenomenological nuclear potentials. The convergence of 
the derivative expansion has been examined by comparing two potentials generated at different 
energies. We have found that the discrepancy is small in the low energy region < ^cm ^ 45 
MeV, which indicates that the derivative expansion works. For quantitative applications to nuclear 
physics, 2+1 flavor QCD should be used to generate the nuclear potentials in the light quark mass 
region. By using PACS-CS gauge configurations, we have attempted to obtain 2+1 flavor QCD 
results of the nuclear potentials. Qualitative features remains the same except for the enhancements 
of the repulsive core, the range of the attraction of the central potentials, and the strength of the 
tensor potential. These nuclear potentials have been used to calculate the phase shifts, which behave 
reasonably. Although they are reasonable in a qualitative sense, their strength is not satisfactory at 
all, which suggests the importance of the lattice QCD calculation in the light quark mass region. 
Finally, we have applied our method to the hyperon potentials such as N'E and A^A, for which only 
a limited number of experimental information is available for the moment. 

It is interesting to use our nuclear potentials to study the nuclear many body problems, which 
provides a way to access nuclei based on QCD. Needless to say, there is another direct way to 
access nuclei by lattice QCD, i.e., direct calculations of nuclear spectrum, matrix elements, etc 



[21]. These two approaches are considered to be complementary. The former keeps a connection 
to the conventional nuclear theory, while a number of formalisms and techniques have to be estab- 
lished. The latter loses a connection to the conventional nuclear theory, while many of the existing 
techniques in lattice QCD can be used. It is desirable to use both of these two approaches as the 
situation demands. All that is certain is that the lattice QCD will provide a unique tool to study 
realistic nuclei in the quite near future. 
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